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ABSTRACT 


The distribution function at a fixed’ time of a random walk 
discrete Markov process is compared to that of an associated 
continuous process. The difference is shown to be proportional 


to the square root of the time increment of the discrete process. 
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CHAPTER ONE 


ESTENSION OF THE BERRY-ESSEEN ESTIMATE 


TO CONVERGING SEQUENCES OF ABSORBED RANDOM WALKS 


The intent of this thesis is to add to the understanding of 
the rate of convergence of discrete Markov processes to a continuous 
process. The main result addresses the fixed time distribution 
functions of convergent discrete Moran type processes. 

The first section of this chapter gives definitions and 
background for some questions of convergence. The second section 
introduces the "boundary value problem" approach to problems of 
convergence. In the third section some previous results are cited. 
Finally the fourth section states the main result, which will be 


proved in Chapters Two and Three. 


Section One 

In this section Markovian transition functions are defined. The 
results Of the thesis are all stated and proved within this context. 
Following the definition the connection with a Markov process regarded 
as a stochastic process is explained. 

The section concludes with some examples illustrating convergence 


problems associated with a sequence of random walks. 


Definition 


Let x be a measurable subset of R, and S a subset of 
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[0,°) which contains 0O and is closed under addition. Denote by B 
the Borel G=algebra of X:. A function { P(x ult Enis pk cee ee Bt 


1s Called a transition function on the state space xX i£: 
a) pore) defines a probability measure over B. 
2) P(x, {x}) =e ae 
Se yr, U) ets = &) Measurable, tor fixed wi eb. 


a) P A\*,°)” satisties the Markov’ identity: 


Pee) i) Pi(y,U) P(x,dy). 


It is possible to construct a family {ze lectern tes sO} Ok wrandom 
Variables on the state space xX such that for each xe xX, (ze nc =O} 
x 


is a (Markov) stochastic process whose starting random variable Za a3 


almost surely the constant random variable x. We refer to this as 
"the process which starts at x". These processes are Markov Processes 
in the classical sense (ci. 5S “Ch. 3), Their joint distributions are 


uniquely determined by their distributions, given the Markov property. 


These distributions are given by the transition function as follows: 
x 
Prob (Zi € U) = Pi (x,U), iB ee Os ees 


Before further reduction in generality is considered, two 
classes (discrete and continuous) of transition functions will be 
Tdentitieda. A transition function Bay a) will be called discrete 
if Sand =X ‘are discrete. A transition function P “will be 
called continuous if S = [0,") , xX is comnected and the following 
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(0) lim = Pe Gcydy )e =205-, 
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ly-x|>6 
: : al 
(i) em = (y=x) PB. (x,dy) = b(x) A 
1 2 10) 5 
|y-x| <6 
ab ; 1 2 
(( ai st) lim = (y-x) P(x,dy) = 2a(x) 
= 0) | y-x| <6 


The quantities b(x) and 2a(x) may be interpreted as “infinitesimal 
mean and variance". Kolmogoroff observed that, if a continuous 
transition function P (*,°) has a continuous second derivative at 


x € X, then it satisfies the "backward equation", 
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B) ) p) 
—- P (x,U) = ——P (x,U) + =P (x, 
ae , U) a(x) we , U) b (x) '% U) 


mene 1, = © 
(en lll pee 299) 
Situations occur where a notion of convergence of discrete 
EGanclelom fLUnCELONS tO) a COneanous) transition tunction 1s. usefull. 


Perhaps the best known example occurs with symmetric random walk, 


described next. Let n be fixed, 
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Note that the last expression together with the Markov identity 
completely defines PY (x, U) for all tes". As n> (for 

those n where x € es and, tré cae application of the central limit 
theorem yields that the distributions of Bean converge to the 


Normal distribution: 


aL ar (vex) °72t 


P, (x,U) dy s 


f ; 2 : 
Note that this choice of T = l/n causes the variance of the 
ela : ; : : , , ; : 
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normalization for the central limit theorem. 
Also the Berry-Esseen theorem gives the rate of 
covergence as L/nvt fetes (2) 7pec0O). teers this mete that 
one expects to find in related covergence problems. 
Observe that, even though "infinitesimal moments" do not exist 


7 n " . . ° . e ja 
Bor discrece {Pi *)); approximate infinitesimal moments" exist as 


indicated below: 


2a_(x) = = ) (y-x) Beixs ty!) = 1, 


n _— 
b_ (x) ) (y-x) 2 1Gy bere” 


In the above expressions, integration over a discrete measure becomes 
a sum. Now it comes as no surprise that the limit transition 


function P (*,*) satisfies the backward equation: 
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Now another example, the Wright-Fisher process from population 


genetics. 
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As before, the one step probability defines Daa) EOr Lai 7 exe oa 
The combinational nature of P| (,0) makes it difficult to analyze 


for large n using the Markov identity. However, the “approximate 


infinitesimal moments"! can be calculated: 


al 2a 
2? = — = = a3 
Us) ies 1 Sale pan (Ca ulgy) ees ro oi 
b (xy) = =F) (y-x) P-(x,fy}) = 0 . 
n T ca 
Then if the right kind of convergence applies Pee) may be 


analyzed by considering solutions to the backward equation 
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In this case the time increment T has no meaning in the 
model giving rise to Po(*,*). Here T was chosen so that 
“approximate infinitesimal moments" converge. 
Two types of situation where this " 
serve mention here. 


The above examples are of the first type. Suppose a discrete 
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Evansitetion function is defineduby tts first step transition 
probabilities. Attempts at analyzing its behavior after many steps 
by using the Markov identity may fail, due to the complexity of 
this one step transition probability. However it may be well 
behaved in the sense that, the transition function may be 
associated with an element of a sequence where the “approximate 
infinitesimal moments" converge. Thus analysis of the original 
transition function may be carried out by analyzing the associated 
backward equation. 

Another type of situation is the reverse of that previously 
described. Suppose a numerical solution to a backward equation is 
desired. One possibility is to invent a sequence of discrete 
Eransition functions, with the first; step transition Lunction! being 
simple (as in random walk), and the “approximate infinitesimal moments" 
converging to the respective coefficients of the backward equation. 
If the right kind of convergence exists, the solution to the 
backward equation may be approximated by choosing a discrete 
transition function far along in the sequence, and then using the 
Markoy identity. 

It may happen that the sequence of processes converges in 
a reasonable sense, and yet the limiting process is not the unique 
diffusion with the limiting infinitesimal moments. Consider this 


example. 
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Let P_(x,{y}) = 1 - Wat fom sya = 0X 
@) otherwise 


Observe that 2a_ (x) Ss it éigel Jo Gx) = @, 
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The limiting process is a spatially homogeneous jump process, and 
is certainly not the standard Wiener process. Note that property 


(olor page, 2-fadils for this process. 


It is natural to ask the following question. Suppose, a 
continuous transition function Ea) satisfies a backward 
equation. Suppose, a sequence of discrete transition functions 


n 4 : , n 
P (*,*) is given. When, in what sense, and how fast does ee rans, 


converge to eee ey, 2 In this thesis, a special case of this 


question is examined. 


Section Two 

beter (>, 1) be a eransi tion funetronsast defined ai he 
previous section. Following W. Feller and others, consider for 
a (temporarily fixed) bounded Borel measurable function f the 


expression: 


p (x,t) - fey P_(x,dy) , xe [0,1], t20. 


(This is the expected value of £(x,) in the process which starts 
at 35) INfoNEEY Velaehe SF GF She = (the indicator function of a 
Borel set JU), then (x,t) = P(x,U) -) Le the “difttusion 
conditzons. (0), (i) (iil) are satistired it follows by a standard 
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= x7 Oh = ob (a) Ene aa (x) ee (ee), a.e. 
Because Byte) is oS (the point mass at x), it follows from 
(4) that v0)" = f(x). “Lf the endpoints “O° 1 “are assumed to be 
absorbi slat P OFT = , e = 
ing en all ) 85 and Pd, ) 51 Keech ie 


Equivalently, the probability of moving away from O or away from 1 
TSeZeOroOw el otsmrorces  UCOeC) aor tO), seul ce = srt lian alte ot. 
These properties characterize the transition function, as 
noted in the appendix. A precise version of the boundary value 
problem satisfied by W(°,°) is stated shortly. 
A crucial, problem jis) to show: the existence of FP (* 7,2) 
having gaven dritt and variance Velocities b(*), a(*), ‘and 
absorbing ends. This thesis needs to use not only the existence, 
but cangexplicit formula, tors P(°,*), in, the case when lIg(2)) = O, 
and a(*) is Riemann integrable with “Oe ks a(:) = k5= the 
existence in a general setting is a very deep result due to Stroock 
and Varadhan(cf. [13], Chapter 3). The appendix proves existence in 
the particular case above, by solving a boundary value problem using 
eigenfunction transforms. The connection is made with the semigroup 
point of view, and the diffusion properties (0), (i), (ii) are 
proved, to hold for the constructed transition function. 
The BVP (boundary value problem) satisfied by P (x,u) is 
now made precise when y= [0,z]. If w(x,t) is defined by 


w(x,t) = P.(x,U), then: 
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A discrete boundary value problem version of a discrete 
transition function is obtained simply by changing notation. First, 
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An G@arly result of this type 1s due to FP. Lax (ef. [9] ) 
his convergence result verified numerical methods for solving 
boundary value problems. 

More recently, semi-group techniques developed by Ethier and 
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distributions of a discrete and continuous process. The continuous 
process is a diffusion with zero drift and absorbing ends. It would be 
desirable to hold this continuous process fixed while it is compared 

LO a Sequence of discrete processes. However, it simplifies the 
technical details if the infinitesimal variance is taken as a step 
function approximation to a continuous function. The discrete process 
is a random walk process with zero mean and "approximate infinitesimal 


variance" equivalent to that of the continuous process. 
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CHAPTER TWO 
Itias ‘the purpose of this: thesis to prove 
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where K is a constant independent of n. There is nothing to prove if 
z= 1. In the remainder of the thesis assume z € [0,l). 

In section one of this chapter the above integrals to be compared 
will be expressed as solutions to continuous and discrete boundary value 
problems respectively. In section two, it will be shown that with a 
condition, the solutions to these boundary value problems can be repre- 
sented as contour integrals in the complex plane. Finally in Chapter 
Three it will be shown that the condition for the contour integral sol- 
utions are met, and these integrals will be compared to give the final 
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The above differential equation holds where a” (x) is differentiable 
and -t > 0. 
This boundary value problem can be transformed into a zero end 


condition problem as follows. Let 
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analyzed later. 

The next task is to develop an equivalent notion for the discrete 
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This is the discrete boundary value problem which will be analyzed later. 
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The first integral is a special case extension of operational 
calculus to an unbounded spectrum. 

The following solution to the continuous boundary value problem 
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Observe that each eigenvalue is on the negative real axis 
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Finally since both sides of the previous equality are continuous, the 
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Since the eigenvalues of ne lie on the negative real axis, 
the contour of integration can be chosen as described below. 
Let part of the contour coincide with the contour chosen in the 
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continuous case, ie. A = (l-r ) - i2r, --~ < r < ™. ‘Two complex 
numbers a and do will be found in Chapter Three which lie to the 


left of all the eigenvalues of baa and then joined to close the 


contour around the spectrum, as sketched below. 
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CHAPTER THREE 


Pevis= the purpose or this chapter to show that 


| 6 (x,t) - (x,t) | < K/n where, 
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o(x,t) = and +. R, £(x) dA 


and 


> (x,t) = = Spann *ete oar. 


It was shown in Chapter Two that this is equivalent to the main 
result, pending the exact description of the closed contour of the 
Second integral. “Iesassumption 2 ¢€ [0,1) as Still kept. Also, the 
restriction and definitions on page 11 still hold. 

To achieve the above inequality, the contours of both integrals 
will be truncated. Then the integrands will be shown to be small on 
the left portions of the contours, and close to each other on the 
remainder. 


The appropriate bounds will be developed in sections one through 


: ; n : : ; 
five. In section one a and Ry will be examined in terms of 
Green's functions in a general setting. In section two details of the 


Green's functions are developed for this specific application. Section 
three has calculations which give a bound for the left most eigenvalue 
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of A , and hence the description of the closed contour can be completed. 
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Section four contains a list of technical lemmas and their proofs, 
which are later used for the comparison of terms. Section five finally 


has the proof of the inequality which begins this chapter. 


Section One 


Observe that R, has the following representation [see appendix] 
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Section Two 

In this section further information on YQ Arx) and ¥, (A,X) 
can be obtained from the restriction that Ae) is piecewise constant. 
The expression of YQ Ar) on the sub-interval [0,x,] can be used 
to find the expression for YQ (Arx) on the sub-interval [x,/x.]- In 
fact the expression of YQ Ar) can be found on any sub-interval by 


finite recursion from the left. A similar recurrence relation will be 
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found for Yo (Ae) « Comparison of these recurrence relations will 
give a comparison of YQ (Ar) and Yo (AX) which in turn aids the 


comparison of R, and Ze 
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An induction argument then shows h < |b. | 
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Lemmas to be used later are listed as follows. 
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Let Ba a Ne elOr cal Vera 
Suppose G has a continuous Jacobian in its second variable, 
with F, G and Jacobian of G bounded with respect to both variables. 
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In the above expression JG is the Jacobian of G in its 
second variable. The sup norms are in both variables. m is defined 


by 
m= nax4| Ger Neen 


Lemma 11 The finite sequences ao; and 8; are uniformly bounded in 


a | 
2 
ST Cae Ae fOr Ly ene yey = plier 
: Vv 3) j 
Lemma 12 If A= (1-r7) = 92, wath hat < Re A and a = jh = aah 
pce ae n 


t< 4 <n, then the following inequalities hold for some constants 


Koy Ke and n_ large enough 
; n 
i) XQ (Ae* “XQ (Ars) |< Ke|A|/n 
nae n 
ia) [¥ Ore) -¥y Oax,) | < K,|A|/n 


_— 


a nl a 


ied el ror 


crint Sy 7 mee: adj eb — @obaerinies atoms ede 
Boar a ks ob "aR rite SSA turn eae, SL “eed eae 


Vy at » Kt + i 
sau 7 sree U 1 


er ee 
2 4-s Qoe feb | (Nigra 4 i" tet wee 
. wa 
; Ay 3 : 
_- =a ® Ce ee eee Se See 
aie iieke te’ 9 ol ets /eurame qainsic® eal madts : 2» 


47. 


Lemma 13 
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Use of the original second order differential equations yield 
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Now that the alternate formula for Rif has been proved, 
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Proof of Lemma 10: This is just the vector form of lemma 9. The 
reasoning in lemma 9 and 10 needs to be adjusted to suit the vector 
forms of the statements. 
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The product formula for checking the vector form of lemma 9 
becomes 


where ee Ls) a matrix acting on on 
The manipulations used to prove lemma 9 applies to the vector 
case if 1 is replaced by the identity matrix. 


Thus lemma 10 is proved. 


: 


_ 7 
of 


on - 

up ws ha 
eee or “ | 

& | ee ore | 


ot! ata spree 
posoay (afd tats odo 


hy i” (Ae: en 
a a 1 


. | - ; 
no paiaeiego ota 2 AE ie hy! a 
Diath I oll 


“ sornev Geely te qth gregs 


tsi 
nosey Ad 10 eatgon pat eet 4 caer 
| pe a by if 


@ nea 30. He tte =" » (leeks 39 ae rat _ 
a, 
‘ a 
ole + 
pee + cepP gee! = (yale ; 


~— 


> 

es eahawe sistem 6 a 

jedoay sit cs apilegs & eum k avouy. oo bear garods 
Estim yo hghebe way, Viana at: 


foe a 


6S. 


Proof of Lemma 11: Recall 


e 
I 


‘ —X 
Fader tee) idk AL Srey es 
‘ 0 | =.) sin /* AEE. Ge es Sa) cos 2, a 
Bee VG ee) eos a (A, < ) sin /— x : 
J 0 J dX ax Yo! a. 5 


MILO Wo S tw sb al, 9 K se — incl © < a < He Furthermore VAY (AX) 
da 

and ax 19 'Ar*) are bounded (from lemma 2). Thus as and 8. are 

bounded. if 


Oh 
Pigeoe Che iene Ws IS IO Wanlil los wisecl, ilbeig 2 ye( ’) be the 


recursion relation taking 


i cheat 
n n n 
o = | 5 Oh 41 
Let F i be the recursion relation taking A TO = : 
Nn A 


The expression of F at the end of section 2 is not immediately seen 

to be close to F. However, if the correct terms are added and 

subtracted to change x, + to se and then formulas for the cos 
J 


and sin of a difference are applied, then F" takes the following 


form 


(ie Adee wiokieadina pee cont ey i fre 
a Le Bust? prs aicial a heeelial wean 


, @ 


gid 6d «a i) +01 tna: ck Thee, OF spite 
pitenay saabeses 


6 


- n 
AN Ge Mad f j.° “_ 
. AM ) oh - ; parva (ox elpr moieties Stops (ded 
# Waal 


4 LF; 
ior Slhternemmsr2cd Al & Hoiosm 1o bro Gd SA 8a Yq SGreevineh is 
tr Dobbs S18 Smead towiagy sis it ad <2’ oF asdte 4 
ét> eis et. eotupd? Hands Bde ve aa £4 a aes. 3. 
seaseptitet eis éox07 9% munl) ~tettggs wae oddities! 50) aata 7 
“” ee 


*\ Goan 


= ’ 


6G. 


n n n 
on = || On Gosh, ith, sai ess scala, lib, - Sib 
jae 5 a ies j 1 gyri Af 
2), gard 
- J -l1](a. cos Sie sheicks Gali Ws GOS Wy, . 3%, 
ok 5 ee mee) pe) sae? a 


n Te n 
on = |) @ Sasa, il, .. owl, ets, Cos G<. lu. ..=, 
gal J a 2 me iam ae 
a n 
+) 1- es (G, Gos aioe ea Sug Ww) oe ,) Saki WM, . Fu 
a4 J aie shes ee jl 3 


ne n 
+(a. Sin U.x.+o. cos eaeEN COSMU 
i Ee aa ape) Jel 
a D 
1+ A 
-/ ee rh eas = PA (o” cos Tee Sin aes;l) 1+ Z F 
Be alcel | J a) a oe ay =. 


In the above expression EN is the difference operator defined 


by 


2 ‘ 
It will be shown that ||F-F||_ < ch fon (some cy 107 e Sarst 


; : n 
note that the first two lines in the last expression of | and 


n 


Pad 


are nearly equal to the expressions for One) and eae The 


ae i *. Wie pee se a 


bp e Cee: 


_ a 
MEE ee aes dle 


et ii al ill am 3 J a 


- 


Pek Sit wre re nti 


. a ob \ Chee aie eo at Prete tis J Mk gett » ole a 


hentia sod oxo aniem2Tlh =o el A hod 2.9 oho sieht : 


“iQ > expe = bey A a! 
. a - » 


7 ay 


sexe? 40.4 5 pie Her (iota lane 
pits 4 ne a ch paalae ca 


sim yyy Bre ge 5 lenin i ag 


OF. 


: : h 
GSIMEISIV Scala Che la Ga in the denominator to the difference of 


3) 
n ; h A : 
F and F is of order h2i if ve is bounded away from ll. 
j 


However, 
|| < Bn a Wy dee : 


SO, 


fon “h small enough. 


—)\ 
Next at willl be shown that aid =| Le SuOi Omer reo Ml c 


j 


2. Semiseh! faa : g 
Recall u, = h Sin Cae Taylor series expansion gives 
: j 
2 3 
= ely. 
You Sin ey | < 


fon 6) 0 and | ey| my aelbics 


h A : 
and we have already seen that bec is bounded below 1. 


TY say ci eh P 


t 
vn asta ts. 40 FEMS a ee ee ee 
N| Tal. 7 
avo ae laiagte 4a.18e ete = ¢ ne al 

tens ; 

pee: ere eo” aie Ay 

igs 

~ tel en2) 


aot ot > |ys| St Me 2 Oe 


Paes | al 
: ‘eft Pal — yt ST 
Pei - 
- 


68. 


So 


A >| < ch*y Ml pe for some )c =) OL So far 
J 


< ch’ [al (fa]+/8]) 


for some “= 0, if the first two lanes of FO are compared to F. 
The last two lines of the expressions for on and BP will 
ee ane! 
be seen to be small of order ee Uo leee les oe as shown by this 


sample calculation. 


Consider the last two lines of the expression for q 


j+1° 
n 
(Gyn “Sali WO ae COS! WSs. )) Sal Wl, . ad 
ai Bs ap) eS al 
eg) n 
+ t= cos BL oy ery Sec (Gr, COS Wok 8, Swi WK, )) 
ar yea 7 a 
=a) 
n Bag h 
- f : i . aye 
(a. Sin Carte COS Bos, 2 Sala a Cos a 4 Zh) 


The above expression is the difference of two nearly equal terms. 


Note that 


| A 


chlu,| for some cy > 0, 


; h 
[ee Ux] = Sin u.(x.+=)| 
ne! he a 


h 
|cos Eee = Cos us (x.+7) | Senn os for some Ope OF 


‘ae ee 

‘J o3 nessa ‘ate, Fk Sp noha cdl Git 3 
ia 4, 8 Baa ra ah a eean oan 
eld qo mire Re (rate al Tera 


Me, 8 aa haa sank oe 


“ht 
Cae ae Ave eb y ity 
( Pat ate mers an ‘ai ene 


re ja es eds 105 feng nie beh 7 


is iva -) . 
Let ad : a ’ it ~ cH bk Re My iy 
~ itt. | et A Bye ee Oe 2 al ied re ee Beg ha: * ; 


«hatest Le1po vizsen- ced Yo wtharettth eAtgL NosaenAy aime 
on ele 


2953 omoe. x07 mes ibaa ee 


as 


= t 
i wes z 
- Ce, wae) 

=~ 


me 


Ga 


and 


2 
Thus the terms being considered aressmalleot orders hey) (locl+les |). 
Recall that lemma 12 says that (ja.|+/8.1) is bounded for 


iv Sve SO 


This gives 


(j-1) ||F’-F|| < he, ||- 
Now consider the Jacobian of F(x,*). Observe that F(x,°) is 
On 
a linear operator on rs and its matrix form is given below. 

3 

cos — v-X x. A os =sin = y=’ x. A == 

a, 
J J 
sin - v-A x. A = cos, = y= x, A = 


0 we 

vitairfinty TANS: vsaeacta ‘ftesit es taailletc mat 
ui Bete el A], abel el) Sete eyng SL, same asd sts 7 
iy A oe mee 


HAL aad = ||) oat 


TY (+299 Sard svasedto «4 »%)% ta ftatddoel dad aber Ge - 
Woted. tevke 28 maio8 Kradem ett. bie. {sd ne conacnap sonal 3 


705 


Ava. Waele J 3 
+ ij 
Va, 
J . -r -r f -A 
O Sin - Ss cos a, -sin Aes 
Bit J J 


The first term in the above expression is a rotation matrix with 


; 1 . ; 
angle dominated by A as Also recall that the imaginary part of v-A 
5 
148) loowbuaclkecl, “UNaWIS| Tela sealteSie (eee) ale) Ib 4p De. where le, is bounded, 


and I is the identity matrix. Also the second term of F can be 
expressed as ae with SS bounded. So JF = (I+hE) with ||| 


bounded. Using the notation of lemma 10, observe 


eee < (1-h|/B|])” 


and 


jae f2> < astnlfel[?. 


a le 
S2nce (= |e i) "and (1+|[=||n)” converge to elEll they are 
bounded. Now consider 
fe} an 
ee a 
B n i 
1 By 
u 
Gy i a 
Ohya ee = ee a . 
= 2 
1 r a8 ihe 
4a 


e So ee 


. oe wad Tl, ns votaviion "Gulaiieee bos "etal 


Patna ei 13) aes ot * 1 atm ae 
ager 1 Do, wie seal ne atk el som 
jaf aéw iis = 36 oe ‘a 7 


wrzeede OD ans a 


om et it)» a 
SAS 


J pai > FP8 jan) 


—o 


ssbluios wal ; ale 


TAN 


2 
Thus Jo,-ai | “sch A forssome. "cr 0. 


n 
ANAL = = (0. 
so By By @) So 
a a” 
1 
- - < BREN LOGNSOMe eG On 
B n 
ii BF 
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lemma 12 (ii). 
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Section Five 
In this section the proof of the main result on page 11 will be 
completed. As previously noted there is nothing to prove if z=l. 
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APPENDIX 


This appendix contains results of two types. The first type 
of theorem would be standard if some restrictive conditions were 
imposed. An outline of proof is given which hurries over the familiar 
arguments, 

The second type of result found here is a direct consequence 
of some very difficult and lengthy work by others (for example, 
pointwise convergence theory, and harmonic analysis of the sort due 
to Elias Stein). Such results are proved here in a form sufficiently 
reduced in generality that elementary arguments apply. 

The first section deals with the spectral theory of Sturm- 
Liouville operators. The second section solves boundary value 
problems using the first section, and the third section proves the 
existence of certain diffusions, used in section 2. The semi-group 


point of view is also briefly discussed in section 2. 
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§1. A class of Sturm-Liouville operators 


Let a(*) be a Borel-measurable function satisfying 
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Proof is standard, using the contraction mapping theorem in C[a,8] 


to obtain local existence first. 


Spectrum, Resolvent Operators, and Greens Functions 


Write w(dA) = ¥5 +1) and note that w(A) = 0 iff A as an 


eigenvalue of A_, with eigenfunction Yg Are). Since A Seal 


2 2 


symmetric operator its eigenvalues must be real, and its eigenspaces 
mutually orthogonal. Separability of ie forces the set of 
eigenvalues to be countable, and obviously each eigenspace is one- 
dimensional. Clearly O is not an eigenvalue, and yet the eigenvalues 
must be non-positive because (a,£,£) Ojala te D,- (Proof uses 
integration by parts.) Although the following lemma is not needed for 


the proofs which follow, it is used in Chapters 2,3, and simplifies 


our presentation now. 


2 Lemma (ct. [8S 1.401 ) 


The eigenvalues of A form a decreasing sequence 
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{A;[J=1,2,.-.3 c (-~,0), satisfying [5 | >cH for some © 2.0. 
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so 


— ° ==: e 
lo) = Ix, 5, ns Yop (Ayet) 


For A ¢ WD eae ote define 


F Uy » ® as 
U ' : DE ad is 20S aX. “ a 
a ; 


: “my, be 4 ‘Ae . 
ne et ( 2y/Or (OR si m4 & ia . 
| ab R saints set) & aoton ee Ae i rm noel 

ge=wyeivpie att Ore hget od sh osigastis ea aonexege Ob at 
a Siig) ond parr 2 2a ‘elhieneet Ssameseose 


-ery 2] seqanarghs faa sae. hae pen ances ad of 7 


seu tavatgie sfid day bra aolarncols ne, Some O yltnaei9 ,tenosm 
Qa fie «02 (? hfe Actin sole Saeinaae ot ‘eumt 


"Ss 
i 
yo? Pabeon tor et snmad print a ens outet La {.ettsa va ue fi 


es 
rey 
a 


a 


eon Poors) * 


fi8 exeracdd n) Beer et ot walint Hobdy atoosy 9 
we ne 

J a = 
(fom.g Led Sev atime tek: 


Sat pes Prekaseow eG MIT: . oA Yo genieynenie att 


as ; . — 
0 <5 oop ~6t C75 } AI wiiviotoeg «(vd 2 benevSsdet} A) 
Ca 


sor \itatia oe 


4jse Bereigorrs Weiaminwle toaibention uit seh eat 


oa 


ip 


seg Mgt” glia gyAhgeR = Mahe 


a he 


83". 


a(y) Ww) ¥) Oux)¥, Oy), x < Y; 
G, (x,y) = 
Si = 
aly) “w(\) Yo Qryv)¥, Asx), xy. 


Define the resolvent operator R,: in a i, by: 
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1.3 Lemma The operator Ry maps Le onto Do, and equals CE 
Proof: (Cutline) Writing out = (R.£) (x) as. the sum of two integrals, 
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one may differentiate directly to show Rf € Do, and (A-A,)R£ = 42, 
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Proof: (Outline) It is easy and standard to show R, is compact. 


r 


Also, one can see that the adjoint of R, is associated with Co a 


and this leads to the conclusion that R, is a normal operator. The 


spectral theorem for a compact normal operator (cf. [14], p- 345 ) 


now yields the theorem. 
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§2. Solution of boundary-value problems using eigenfunction expansions 


2.1 Theorem Let a(*) be measurable on [0,1], with kK, ENGR) Mee ale 


2 
Let fe Llo,1). Then there exists a function ¢(*,*) continuous on 
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therser, [O71] x (C,°)) satistying: 
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3) oo Lsecontanuous on | LO pl x (0, <<). 


3 
4) = Ls absolutely continuous on {0,4}. Se >) 0; 
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6) If £€ [0,1], then lim, ||¢(+,t)-£(-)||,, > 0. 
0) 
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the series converging uniformly on [0,1], for t > 0, where £ =o 
J J 


Note that 7) includes uniqueness of $(*,°%). 

Proof: Consider the series of 7). It was proved in Chapter 3 that 
ale 1s bounded by a power Of J. Lt as also a standard result 
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Since a OU Ge) = Me Ke) asses of [0,1], the series 


eo et 
a(t) ) a Y- (+). converges) in the norm of L [0,1], and it then follows 
that the series is almost everywhere equal to a/(-) “ jmeare On eo neat ace 
ox 


that the solution must be given by 7) follows by integrating 5) against ¥.. 
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Properties 2), 6) remain to be proved. For 2) note that: 
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Property 6) is considerably more subtle, as the series for f need 

not converge uniformly. Its proof will be included in the next section 


on semi-groups. 
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Let £ € C,l0,1], and write Tf for, the solution —o(*,t) of 
the BVP of the previous theorem, when t > 0. Write To = I, and note 
that uniqueness of the solution to such a BVP implies that 
= T T . We will prove that Cr is a strongly continuous 
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the previous theorem. In order to construct a Markov Process from bey 


we need properties 5), 6) below. Property 5) is obvious, and 
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broperty 6) will be proved with. 3). 
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This iS an easy consequence of continuity of 9» on [0,1] x (0,%). We 
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— re) 
D, sup (x,t) = sup (52 (x,t) |x}, 
x 
where Ge is the set where ¢ assumes its maximum, and D, denotes 


the upper derivate. An analogous argument shows that Dy sup (x,t) 
= 


no [1,") to tends Sosale sie so: 4 er oe 


Pr} ($,%)4 ‘Sao oer aaa eamuees 4, T fo briee 


= tl j 4 (2 
. 7 io 
io <@ xed ted9 mtalo 
- _ 
— 
n 
of : 7 a0 ay * a > r * 4 > 7 
‘on : Fax) (0%) iia = (Sox bat 42 
x 
i 
ay vr ri. | ne @ 2 wtunieieg- To andetresios YHso Ae % 


: a _— 
Ayitmoen taces 4& Wo ty Ede VE oePoreer tans wt fe avery won 


iy det oO > {SyxbP Re ee | oo Biel Tet Uo avig amoe 38. | 
, . | 
5 we 1 nol 
ewer F >. .% ihe tet ovizecet 6a JeOW 125.) SZ awoda 
“ v 4% ~ 7 
q 7 > 4 " 
eee 2 ry = (2%) Fo (hie Bumte «,, *) Se [erteial ms eo. aviter a . 
a oF eh ’ Nd a 
= 
jumps, LGvtes0: mB av av svr2spen a2 (ts) ae ebcrbanon ov 408 


: Rs 4 _ 
f2tcr4 ANtIoidewInes SAT . 42 .* nolwplioeG Shy Ret nica abs Be 


(a 21 doit jsviosgem se vnr ek we 
GO * 2 to talo Woe cen ov  tfanraps py sae ankle ee 
a ee 


re 
2 


7 jaail tour) a) we ae fi 1 


> 


FF 


as ait sen 3 ah 


ih 

eh 
a - - “ : = 
A a ’ _ : 


> 
> = 


—— 
‘, a 9. 


ie 


</> 
. 


88. 


must always be. < 0, and hence sup 9(x,t) is non-increasing with time. 
x 


This proves 3). 


Proof of 4): Because the operator norms ea t> 0 “are uniformly 
bounded, we can show that the functions f for which 4) holds form 


[O,11. Tt is easy to show that 1f — 2s one 
Kate 


of the eigenfunctions wa then 4) holds, because ee =e J Hoe 


a closed subspace of Co 


It remains only to show that the closed linear space generated by the 


eigenfunctions is all of C,l0,1]. Tor prove this, wise the Green's 

function representation of R,, to show that [heel ees K|| £||,- This 
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means Ry is continuous as an operator from Li to C,l0,1]. Since 


: ‘ ! 2 f 
the linear span of the eigenfunctions is dense in Li and is mapped 
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dense in the range of R,- But we already showed that the range of 

Ry is D., , so the final step is to show dD, is uniformly dense in 
Cyl0,1], which is a consequence of Lemma 3.2. 0 
Extension of Ty ee ClO, 1) 
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Cj lent] : 


ani puna olhge th wrron aotaroae & 


ore? ablori : (asta 1 > inate | fala 


So 242i (>. 5% gti Made 99” “ena mt ot NEON ES 
# i 
.¥ eae ets? marae ebins (+ here 7 Sal i minietaae i be 


“sily va bs 1p aanen ones apoult baaok: a eb oe whey ese 2 
= nweto SAT ony El Ap ievesg, of T1073 ey bis al en 


eit: +. 43 1m ya fil dmitd wort ot 4 FO neigesaeersge7 6m. 


oomt2 [E019 of *2 mor? is slated be va a6 ewatmisnon ai - 
Letina. &F- Gop a Ai paneh el ere rerre ens lo ne{e raonkt @ 
\ daotiaw at enoijoabengpie SAs 5p Cage TEATS aiid “ff Maes oF | : 
lo wblres of) secy bowie: YOnet he. ov te . pe io apres off nd eure, 
7 


rk) oem yyEgo* Ta G2. <M Weds os "ad gore anit est om, oO ef ' 


f.£ nmin So-eyneigecnos se at dpiav {£0149 | 
e v1 >» q 


_ 


(£yO3R9 mat a” to 
Fels 


——— 
He 


ot OL) bokkv@ye ebind ifase 1.8 mexoaat pads y [2010 @ 2 iz . aa 


ebiveup oF - GS ® . IR = FL) h , (OD = foe 1D we boca = 


~ 


% ‘Bee J toolrouut amgee owt std 10% daldoxe add aviow any 


envien 049 [8 Ge = (a\e)4 eto dade | oT tot 1 


tol drut 2) tone bedewe (bs -, 0S inieehe niente 


89 .. 


§3. Existence of diffusions 


Suppose 0O < kK, Sk Bd Sd. as before, and make the 


2 
additional assumption that the set N of discontinuities of a(-°) 


has measure zero (i.e. a(*) is Riemann integrable). 


3.1 Theorem There exists a transition function P{-.*)) shaving 
the properties: 
: -1 : A 
(o) tim 1c P(x, [0,x-5] Uc Isto ly) =e forcakiy hy 0p ce Ont ie 
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: sl 
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(iv) P (*,-*) generates a strongly continuous semi-group acting on 


Elor ll. by tne formula: 


(T,£) (x) = f £(y)P, (x,dy). 


The rest of the section will outline the proof of Theorem 3.1. 
Theorem 3.1 has essentially been proved (in a different form) ina 
very general context by Stroock and Varadhan (cf. [13], Chapter 3), but 
the proof is extremely long and hard. Because the spectral theory just 
developed is necessary for the text, the opportunity is taken to use it 
to give an “elementary” proof of Theorem 3.1. 
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and the associated BVP of $82. The solution $¢(°*,*) which was 
produced gave rise to a semi-group {T, | t20} one Cl0,1], as we 
showed, where (Tf) (x) = Ok, tee rae trea lle ee ei oe, ahie 
observe that the functional f > (TL £) (x) is positive with norm 
bounded by 1. Since T 1 = 1, the norm of the functional is exactly one; 
It follows from the Riesz representation theorem that f is 
represented by a unique probability measure denoted Pi(x,*), in the 


sense that 


Cr) =) hc) ee Gcey)” 


We assert that Py (x,*) restricted to (0,1) is an absolutely 


continuous measure. To see this use Theorem 2.] to show that 


eee 3 él, for £ ¢ C,[0,1]. 


If the measure were singular, Ty could not even be bounded. 
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Proof: It suffices to show that there exists for any @ <8 a 


function £ satisfying (i), (ii), together with: 


(W)C) Oe El) 


(vy —£"(B) = 0) £(8) "= 2. 


Choose a continuous function “~h” supported on ja+e,;6-e]” satisfying: 
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h(x) < 0 a S (at8) : 


A -1 
f hlx)a(x) “dx = 0. 
@) 
x 
Write g(x) = if haya (a), “du, 
@) 
sie 
f(x) = { g(u)du 
0) 
TE 3S Gasily seen that a(-)i"(-) = ht) a.e., and all the remaining 


properties can be made to hold by multiplying f by some constant. 0 


Incidentally, this is the lemma needed to prove D. is uniformly 


dense C, [0,1] , hence also dense in ibe 


3.3 Lemma Suppose a(°)f"(-) is equivalent to some function 
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Lins {f £ly)P, (x,dy)-£(x)} = h(x), 
t+0 


unitormly cor x ¢ [0,1]. 
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PFOOFs Write. 6 (x,t) f £(y)P, (x,dy), 


W(x,t) = f h(y)P, (x,dy), t > 0. 

By definition of P, (x,dy) these are solutions to BVPs with 
Starting functions £, h -respectively. If both £ and h vanish at 
O,1, we may apply Corollary 1.5 to the eigenfunction expansions given 
in Theorem 2.1, deducing that Ao$ (+ ,t) SU ec) at On OC aE y ah 
do not vanish at 0, 1, we appeal to the last paragraph of §2 .) 


Then: 
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But both sides are continuous in x, hence equal everywhere, and then 
Theorem 2.1 implies that the right hand side tends uniformly to 
Wx, 0), as = oe Since w(x,0) = h(x), we have finished the proof. 
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3.4 Lemma Condition (oO) of Theorem 3.1] is satistied by ise Gorgias 


Proof: Apply Lemma 3.3 to a function £ satisfying the conditions 
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lim, mu f(y) P, (x, dy) -£(x)} =eav(oe) ee (x0) 
t>0 


at all points x where a(-) is continuous. 
Proof: As in the proof of Lemma 3.3, write 


e ese 
a £(y)P, (x,dy)-£(x)} = t f wlx,s)ds, > 0, 
O 


almost all x. 


This= reduces prook to showing that 2£ a(-) “is continuous: at 9x, then 
lim wis Ss) = alxpet (x). 
s70 


In the notation of Lemma 3.3, this means proving: 


lim i hly) PR ,dy) = h(x). 
+ eZ 
t>0O 
Proof will be finished by outlining a standard "pointwise convergence" 


argument which applies here. Approximate h above and below on 


l’ h,- Use condition (0) (which 


(x-6,x+d) by constant functions h 

we just proved in Lemma 3.4) to ignore integration outside (x-d,x+d). 

Next show that lim sup f h(y)P, (x, dy) may be estimated above using 
to" 

h., and that the lim inf is estimated below using ho, finally 

concluding that the difference is arbitrarily small, hence zero. This 


proves lim f h(y) P, (x, dy) = h(%)), as required. 
tot 0 
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Proof of Theorem 3.1: Properties (iii), (iv) were established in 
52,7) end | (0) “was proved in 3.4. For the case «6 = 1, condition (Ge) 
now follow from Lemma 3.5 by setting f(y) = y, and noting f"(y) = 0. 
A Xo issa POint Of cOntanuity of Jae), write: =(y)= (y-x.)*, so 
a(x )E"(x.) = 2a(xo), and then Lemma 3.5 applies to f£ at Xo 
Vea Gang. sage) ve 

To pass from the case 6 = 1 to general 6 > O, use condition 
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